Experimental study of 113 Cd beta decay using CdZnTe detectors 
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A search for the 4-fold forbidden beta decay of 113 Cd has been performed with CdZnTe semicon- 
ductors. With 0.86 kg ■ d of statistics a half-life for the decay of T 1/2 = (8.2 ± 0.2(stat.)t°'^(sys.)) ■ 
10 15 yrs has been obtained. This is in good agreement with published values. A comparison of the 
spectral shape with the one given on the Table of Isotopes Web-page shows a severe deviation. 

PACS numbers: 



I. INTRODUCTION 

In past decades the investigation of /3-decay played a 
major role in understanding the structure of weak in- 
teractions. Nowadays, such studies seem to be a little 
bit out of fashion, but there are still a lot of interesting 
problems to be investigated. Among them is the mea- 
surement of higher order transitions like 4-fold forbidden 
non-unique decays. Their log ft- values are larger than 
20 and correspondingly half-lives around 10 15 yrs and 
longer have to be measured. Only three isotopes of this 
type are known, 50 V, 113 Cd and 115 In. No beta spectrum 
measurement is reported for 50 V and only two recent half- 
life determinations for 115 In 0,0- 



The transition 



J Cd 



113 In is characterised as a 



1/2+ — » 9/2+ transition with a log ft-value estimated 
to be 23.2 Q (Fig. QJ. A first attempt to measure this 
decay with a small CdTe device resulted in a rather large 
half-life range of (4 — 12) • 10 15 yrs as reported in j4j. 
Two more recent measurements using CdWC>4 exist. One 
measurement using the material as a scintillator resulted 
in a half-life of (7.7 ± 0.3) ■ 10 15 yrs 0. The second 
experiment used CdW04 as a cryogenic bolometer and 
obtained a half-life of (9.0 ± 0.5(stat.) ± l(sys.)) ■ 10 15 
yrs 0. Here, we report on a completely independent 
new experimental approach by using room temperature 
CdZnTe semiconductor detectors. The measurement has 
been performed within the programme of test measure- 
ments for the new COBRA double beta decay experiment 
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II. EXPERIMENTAL SETUP 

The data presented are obtained with four 1 cm 3 
CdZnTe detectors provided by eV-PRODUCTS based on 
coplanar grid technology, i.e. reading out only the elec- 
tron signal. Measurements have been performed in the 
Gran Sasso Underground Laboratory (LNGS) in Italy 
providing a shielding of about 3500 mwe. The four bare 
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FIG. 1: Level scheme of the 113 Cd decay. 



detectors are mounted in a copper brick with all pream- 
plifier electronics moved to about 25 cm distance. This 
copper brick itself is part of a (20 cm) 3 cube made out of 
electropolished copper and a further passive shielding of 
15 cm of lead. The whole setup is located in a Faraday 
cage made from copper plates. The cage is surrounded 
by a neutron shield, consisting of 7 cm thick boron-loaded 
polyethylene plates and additional 20 cm of paraffin wax 
at the bottom. The full paraffin shielding was finished at 
a later stage. 

The energy resolution and stability of the detectors is cal- 
ibrated regularly with the help of 137 Cs, 60 Co and 228 Th 
sources. For the four detectors a CAMAC based data 
acquisition system is used. The signals are fed into four 
13-bit peak sensing ADCs (types LeCroy 3511 and 3512). 



III. DATA EXTRACTION 

The total measuring period consists of 4781 h. Each in- 
dividual run is limited to one hour. However, some data 
cleaning has been necessary for this low energy range. 
Due to noise problems several run periods had to be per- 
formed with rather high discriminator thresholds. As a 
consequence, one detector cannot be used for the low en- 
ergy analysis presented here. Thus 25 % of the total data 
had to be rejected. In addition, some data sets had to 
be rejected because of short term distortions of the de- 
tectors. If the event rate per hour is histogramed, the 
distorted runs can easily be identified by having excep- 
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FIG. 2: Event rate distribution of individual runs for a se- 
lected period of time, shown for the energy range of 100 - 
320 keV. Runs with electronic disturbances would show up 
with rates well beyond 30 events per hour and can be easily 
rejected. 

tional high event rates (Fig. [2J - This has been done for 
every detector independently. In this way further 57 runs 
were discarded. An instability in the threshold leads to 
a loss of 12 days for one of the ADCs. 

For this purpose rather high energetic calibration 
sources were used, thus the datasets have been checked 
independently for background lines in the region of the 
113 Cd decay. It has been found that the 351.9 keV 
gamma line from the 238 TJ decay chain, which is located 
just next to the endpoint of the beta spectrum, gives a 
suitable handle on the quality of the calibration in the 
low energy region. The final set of good data has been 
selected by checking the position of this gamma line for 
each measurement period and taking its position into ac- 
count for the final calibration. The final sample com- 
prises a statistics of 0.86 kg-d. 

IV. RESULTS 

A. Half-life of 113 Cd decay 

The spectral shape of the energy spectrum of electrons 
emitted in a 4-fold forbidden non-unique beta decay can- 
not be predicted by theory, hence we follow the procedure 
described in 5] . A model spectrum is built using the form 

N(E) = F{Z, E)p{E + mc 2 )(Q - E) 2 S{E) (1) 

with F(Z, E) as the Fermi function, Q the endpoint en- 
ergy of the beta spectrum, p the electron momentum and 
E the kinetic energy of the electron. The correction func- 
tion S(E) is assumed to be a polynomial of the form 

S(E) ~ p 6 + CipV + C 2 p 2 q A + C 3 q 6 (2) 

with q as the neutrino momentum. Such a high order 
polynomial provides a reasonable fit and is motivated by 



theoretical arguments on 4-fold forbidden unique beta 
decay. However, the decay of 113 Cd is known to be non- 
unique and it is not obvious that this spectral shape is 
correct. The correction functions used in this analysis 
are based on the measurement of the coefficients given 
by 0. 

The corresponding model spectrum is folded with the 
detection efficiency, calculated with a Monte Carlo simu- 
lation based on GEANT4. Due to the short range of low 
energy electrons in CdZnTe the effect on the spectrum is 
marginal and corresponds to only a modest shift of the 
maximum of 0.2 percent. 

To get a more sophisticated energy resolution for the 
low energy regime, intensive post calibrations have been 
done using 152 Eu, 133 Ba, 57 Co, 241 Am, 60 Co and 137 Cs 
sources. The measurements can be described well by a 
linear dependence of the energy resolution as a function 
of energy in the range of 100 keV - 1.4 MeV. For the col- 
lected data the energy resolution is extrapolated from the 
high energy calibrations done during the measurements 
resulting in an energy resolution (FWHM) averaged over 
the detectors of 

AE(kev) = 4.7% • E(keV) + 28keV (3) 

Before fitting the experimental data with the polyno- 
mial form, a background component is subtracted from 
the single spectra. Therefore, an exponential background 
of the form 

B{E)=B 1 exp{-E/B 2 ) (4) 

is fitted first, well motivated by the observed data 
(Fig. To avoid any contribution of the 113 Cd beta 
spectrum itself the background fit is done using the en- 
ergy range from 400-1000 keV. The noise threshold is not 
included in the fit, it is well below the chosen boundaries 
for the spectral fit described below. Before performing 
the fit, the 351.9 keV gamma line with its Compton spec- 
trum is removed for each detector independently. Again, 
this is based on a GEANT4 simulation of energy depos- 
tion of external 351.9 keV gammas in a CdZnTe detector, 
smeared with the measured energy resolution. 

An absolute calibration of the detectors has been per- 
formed after the measurements. Here, a 137 Cs source has 
been used to determine the relative efficiency of the de- 
tectors. Normalising to the most efficient detector, the 
remaining two had an efficiency of 67 % and 80 % re- 
spectively p|. The major additional uncertainty is the 
active volume. However, various information about the 
detectors provided by the supplier and the observation 
of alpha particles created at the surface lead to a con- 
servative upper limit of regions with reduced signals of 
10%. 

Finally, the three individual spectra are summed tak- 
ing the efficiency of the three detectors into account. The 
fit interval ranges from 120 keV — 310 keV, the endpoint 
of the expected spectrum using the latest number on 
atomic masses is Q = 320 ± 3 keV [j| . The fit range 
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FIG. 3: The summed low energy spectrum of three CZT de- 
tectors. The red curve corresponds to a exponential fit in the 
region 400-1000 keV. 

contains a total of about 37000 events, 12.6 % can be 
ascribed to background. 

Systematic effects have been checked by varying the 
ranges of the fits for the signal and the background re- 
gions, as well as the energy resolution which has been 
folded with the theoretical spectrum. All these variations 
cause an effect on the calculated half-life of less than 1%. 
The slight difference in spectral shape between the two 
best fit spectra of [j| and |(| is negligible, being washed 
out by the energy resolution. As the typical Zn admix- 
ture is between 7-11 %, the uncertainty in the amount of 
113 Cd in the detector material is 1.8 %. 

The combined energy spectrum is shown in Fig.^ Us- 
ing the spectral shape of 0] including their best fit values 
results in a decay rate of (15.98 ±0.41) per hour and thus 
a half-life of 

T 1/2 = (8.2±0.2(stat.)l^(sys.))-10 15 yrs . (5) 

This half-life is in good agreement with the values quoted 
in 0-0 Note that the half-life has not been determined 
by a real fit to the spectral shape, only the normalisation 
of the fixed spectral shape has been changed. The in- 
dividual numbers including only statistical uncertainties 
areT 1/2 = (8.1 ±0.2) • 10 15 yrs, rVndf = 117.2/169 (De- 
tector 1), Ti/2 = (7.5±0.3)-10 15 yrs, x 2 /ndf = 160.4/169 
(Detector 3) and T 1/2 = (9.3 ± 0.4) • 10 15 yrs, x 2 /ndf = 
318.1/169 (Detector 4). 

B. Spectral shape of the 113 Cd beta spectrum 

Trying a real fit to the data in order to determine a 
half-life suffers from a number of problems: Several op- 
tions exist to describe a 4-fold forbidden non-unique beta 
decay spectrum. They are compiled in Fig. EJand are al- 
ready smeared with the energy resolution. The spectra 
of p| and agree reasonably well. Another option is 
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FIG. 4: The summed low energy spectrum of three CZT de- 
tectors after background subtraction. The shoulder of 113 Cd 
is the most prominent feature in that energy range. Also 
shown is the best spectrum of the type given in Q describing 
the data. Note that this is not a fit to the data, only the 
normalisation was varied. 




E/keV 

FIG. 5: Various possible spectra discussed for 113 Cd decay. 
Shown are the spectra of [f| (red), |(| (green), a rescaled 115 In 
spectrum [Iol| (blue), the Table of Isotopes Web-page spec- 
trum flUl (black) and the measured spectrum described in this 
paper (yellow) before applying a cut on the noise threshold. 
The solid curves correspond to the original spectral shapes, 
while the dashed lines include the energy smearing due to the 
energy resolution. 



to compare the spectral shape with the one obtained for 
the 4-fold forbidden non-unique /3-decay of 115 In and to 
rescale it to the Q-value of 113 Cd. In a polynomial 
is introduced to describe the 115 In spectrum. During 
the development of a low-energy solar neutrino detec- 
tor based on an In-loaded scintillator (LENS) the beta 
spectrum of 115 In has been remeasured and their best fit 
values |ToT ] were used to obtain the spectrum shown in 
Fig- El Last, but not least, the spectrum from the Table 
of Isotopes Web-page can be used ^l] . 
The measured sum spectrum obtained with the CdZnTe 
detectors agrees reasonably well with the one expected 
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FIG. 6: Comparison of the observed energy spectrum with 
the spectral shape taken from the Table of Isotopes Web-page 
[TH|. assuming a Q- value of 320 keV. 

from 115 In. However, it would have a significant effect on 
the event number and thus the half-life, if extrapolated 
below 120 keV, because the maximum is below 100 keV. 
Additionally, in this energy region the noise threshold for 
the measurement has to be considered. Approximating 
the threshold by an exponential function and subtract- 
ing it from the total spectrum results in a residual which 
looks more like the spectrum from 5] and Thus, 
it was decided to use their spectral shape, because they 
have already performed measurements of the 113 Cd spec- 
trum below 100 keV. To solve this issue of different spec- 
tral shapes completely, the spectrum has to be measured 
down to at least 50 keV, which will be possible with the 
next step of COBRA, consisting of 64 detectors. 
In addition to the polynomial approach, the electron en- 
ergy spectrum of the Table of Isotopes Web-page has 
been applied |ll| . The efficiency correction and energy 
folding is done analogously to the procedure described 
before. This spectrum approaches zero only in the low- 
est energy bin. As can be seen in Fig. the spectral 



shape is not reproduced in the data. The spectra are 
normalised to the total number of events in the range of 
100-320 keV. 



V. SUMMARY 

A new double beta decay experiment COBRA is 
planned using CdZnTe semiconductor detectors. As a 
side product of the test measurements for this project, 
the 4-fold forbidden non-unique /3-decay of 113 Cd has 
been measured. A first observation using four crystals 
clearly shows the signal, but only three of them could be 
used for a reliable half-life determination. The half-life 
measurements obtained from each of the detectors are 
in good agreement with each other and with published 
values. In the near future 64 CdZnTe detectors will be 
running, allowing a precision measurement and system- 
atic studies of this decay by extracting single independent 
values for each detector. This will also include lower noise 
thresholds allowing the uncertainty which arises from the 
lack of knowledge of the spectral shape to be removed. 
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